Dissipation- Assisted Quantum Information Processing with Trapped Ions 



A. Bermudez/ T. Schaetz,^ and M. B. Plenio^ 

^Institutfur Theoretische Physik, Albert-Einstein Alle 11, Universitdt Ulm, 89069 Ulm, Germany 
^Albert-Ludwigs-Universitdt Freiburg, Physikalisches Institut, Hermann-Herder-Strasse 3, 79104 Freiburg, Germany 

We introduce a scheme to perform dissipation-assisted quantum information processing in ion traps consid- 
ering realistic decoherence rates, for example, due to motional heating. By means of continuous sympathetic 
cooling, we overcome the trap heating, and show that the damped vibrational excitations can still be exploited 
to mediate coherent interactions as well as collective dissipative effects. We describe how to control their rela- 
tive strength experimentally, allowing for protocols of coherent or dissipative generation of entanglement. This 
scheme can be scaled to larger ion registers for coherent or dissipative many-body quantum simulations. 



Among the different platforms for quantum computation 
(QC) [1] and quantum simulations (QS) [2], laser-cooled 
atomic ions in radio-frequency traps stand out as excellent 
small-scale prototypes [3, 4]. In the quest for large-scale QS, 
the development of micro-fabricated surface traps [5] provides 
a route towards large and arbitrary trapped-ion lattices capable 
of simulating a variety of exotic many -body models. In com- 
bination with the ability of transporting ions between different 
zones of the planar microchips, this technology would also al- 
low for scalable QC [6]. Nonetheless, its success depends cru- 
cially on our ability to overcome the impact of various sources 
of decoherence. Of primary importance is the electric field 
noise from the trap electrodes, which leads to the so-called 
anomalous heating [7]. This heating becomes critical for the 
miniaturized traps required for scalability, where the ions lie 
closer to the electrodes and are subjected to stronger fluctuat- 
ing fields, and even their gradients. Therefore, the anomalous 
heating hampers the scalability of phonon-mediated QC [8], 
and spin- and phonon-based QS [9] in surface traps. 

A strategy to minimize this heating is to cool the complete 
setup to cryogenic temperatures [10]. Another possibility is 
to clean the surfaces of the electrodes by laser ablation [11] 
or ion bombardment [12], which has been shown to reduce 
the motional heating rates significantly. Although these ap- 
proaches are very promising, a substantial residual noise still 
exists, and it would be highly desirable to develop alternative 
schemes to mitigate its impact. We hereby propose to over- 
come the anomalous heating by continuously applying sym- 
pathetic laser cooling during the whole QC or QS protocol. 

The technique of sympathetic cooling requires active laser 
cooling on a particular subset of ions in a Coulomb crystal, 
while relying on the Coulomb interaction to passively cool 
the remaining ions. In addition to its use for precision spec- 
troscopy, this technique may overcome the motional heating 
for QC purposes [13]. In fact, sympathetic ground- state cool- 
ing [14] has been implemented between the sequential appli- 
cation of two-qubit gates to recover high fidelities [15]. For 
scalable approaches based on miniaturization, the increasing 
impact of the motional heating reduces the fidelities achieved 
in this way. Therefore, it would be useful to perform sympa- 
thetic cooling continuously during the whole protocol. So far, 
three different schemes have been proposed: ( i) In the absence 
of fluctuating electric-field gradients, one can mediate the in- 
teractions by certain vibrational modes robust to the heating, 
while sympathetically cooling the remaining modes during 




Figure 1. Dissipation-assisted quantum information: (a) Two- 
isotope artificial Coulomb crystal with ions residing in the individual 
minima of a surface trap. The laser-cooled ions (red dots) assist the 
collective coherent/dissipative dynamics of the effective spins of the 
physical ions (blue dots). Both isotopes may be stored in the same 
individual minima without affecting the resultant lattice geometry, 
(b) Proof-of-principle experiment with three ions in a conventional 
rf-trap. The red arrows correspond to a standing wave providing the 
Doppler cooling of the central ion, whereas the blue arrows lead to a 
running wave tuned to the axial red sideband of the outer ions. 



the operation [13]. (ii) By using far-detuned state-dependent 
dipole forces with negligible motional excitations [16], two- 
qubit gates do not rely on the motional coherence, and can 
thus withstand a heating/cooling that is considerably slower 
than the gate speed. ( Hi) For ground- state cooled crystals, ac- 
tive or sympathetic resolved- sideband cooling may provide a 
dissipative force during the gate operation. For gates faster 
than the heating, the success rate and fidelity of these opera- 
tions can be substantially improved [17]. 

We emphasize that none of the above requirements are met 
in state-of-the-art surface traps, (i) Since the ions lie close to 
the surface electrodes, electric-field gradients avoid the isola- 
tion of robust vibrational modes whose excitations can medi- 
ate the interaction, (ii-iii) Heating rates in room- temperature 
setups (of the order of 1 phonon/ms [18]) coincide with the 
speed of far-detuned [16] and dissipation-assisted [17] gates, 
compromising the success of these protocols. In this work, we 
propose a dissipation-assisted protocol based on an always-on 
sympathetic cooling, which provides a different approach to 
overcome the anomalous heating. We identify regimes where 
the sympathetically-cooled vibrational modes can be either 
used as mediators of coherent interactions, or exploited to tai- 
lor the collective dissipation for quantum information process- 
ing. Since we only require to remain in the Doppler-cooling 
regime, this proposal can be applied to larger registers for QS. 

The model- We consider an arbitrary array of A^t =Na^NT: 
trapped ions in equilibrium positions {rf ^rj} (Fig. 1(a)). A 
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fraction of them, A^^, is laser cooled via a dipole- allowed tran- 
sition |g) ^ |e) with frequency COq and decay rate F^, provid- 
ing the sympathetic cooling of the remaining ions, Na- Two 
hyperfine ground- states {It)^!)} form the spins/qubits for 
QS/QC purposes, such that their decay rate is negligible, and 
their energy spacing is COq (ft = I). To minimize the action of 
the cooling laser on the spins [14], one may use beams tightly 
focused on the T- species, or exploit two ion species/isotopes 
with sufficiently-different transition frequencies. 

By tracing out the electromagnetic bath, we obtain the mas- 
ter equation p = — i[if ,p] + ^(p) with the Hamiltonian H = 
H^^Ha +^ph + Va.ph + Vz^ph, and dissipator ^(p). Here, 

I i n 

represent the laser-cooled ions, the spins, and the vibrational 
excitations, respectively. We have introduced = \ei){ei\ — 
\gl){glV of = |t-)(t/| - \ii){Ul and the creation/annihilation 
operators bl^b^ for a particular phonon branch [19] with fre- 
quencies COn. The spin-phonon coupling required for QC/QS 
is provided by a pair of laser beams in a Raman configura- 
tion [20], which drive a red-sideband interaction 

in 

(2) 

where = |t;)(ii|- Here, Q^a > is the two-photon Rabi 
frequency, and C0a{\^a) is the Raman frequency (wavevector). 
The Lamb-Dicke parameter r]^ = k^j • Cp/ y/2m^^ describes 
the laser coupling to the ^-th normal mode with displace- 
ments along Cp, where rua is the ion mass. This spin- 
phonon coupling arises from the dipole laser-ion interaction 
after expanding T]^ <C 1, such that COa ^ COq — COn, and setting 
Q.a <C cOn to neglect the contribution of the carrier terms. Ad- 
ditionally, we laser cool the T- species in the node of a standing 
wave which is red-detuned from the dipole- allowed transition 

^T,ph = X^^",T+e.e-^^^^ +H.C., = ^ri^^in. (3) 

In 

Here, tj^ = \ei){gi\, Qn = b^-\-bl, and the remaining param- 
eters are defined as those below Eq. (2), but applied to this 
particular transition. The differences between Eqs. (2) and (3) 
are due to the different regimes Ffj <C (0^ <C F^, which do not 
allow us to resolve the sidebands of the dipole-allowed transi- 
tion. Note also that the carrier term is exactly cancelled only 
at the node of the standing wave. Finally, the dissipation is 
described by ^(p) = ^o(p) + (p), where 

^O(P) = I^T (Tf PT+ - T+T-p) + H.C, 
/ 

is the usual dissipator for a two-level atom at rest, where we 
have considered that the typical ion distances lead to negligi- 
ble collective effects in the spontaneous emission. Finally, 

^l(P) = I^5r^,^mV(G^pGm-2nGmP)'r;^+H.C., 

Inm 

describes the recoil effects, where we have introduced F^ „^ = 
2T^aqT]lT]l^in^im. and (1 + 3^^)/10(l +^^) depends 
on the linear/circular laser polarization ^ = 0, ±1 [21]. 



In addition to the Doppler-regime condition (O^ <C F^, we 
further impose that <C (5)^,^1 t- In this case, the laser- 
cooled ions reach the steady state very fast, and can be inte- 
grated out from the dynamics [22]. This can be accomplished 
by a Schrieffer- Wolff transformation for open quantum sys- 
tems [23], which leads to the following master equation 

A = -i[ifcT+i:^ph + Vc,,ph,Al] (4) 
where jl = Tr^{p}, and we have introduced the dissipator 

n 

Here, the effective rates are expressed slsT^ = S{^C0n), where 
S{cOn) = Jq dt^Q'^^^' {Fn{t')Fn{0))ss is the steady-state fluctua- 
tion spectrum of the sideband terms Fn = ^iVn-^ini^i^ + 
Tf), evaluated by means of the quantum regression theo- 
rem [22]. At this level, we can introduce the effects of the 
anomalous heating by substituting F+ ^ F+ = S{—(On)-\-r^, 
where F^^ is the anomalous heating rate for each mode. The 
cooling rates and mean phonon numbers are 

Wn = Re{F- -F+}, fin = Re(F+)/W, (5) 

Therefore, we can overcome the anomalous heating by shap- 
ing the fluctuation spectrum of the laser cooling force Wn>0. 

The master equation (4) describes an array of spins coupled 
to a set of damped vibrational modes. The idea now is to 
use the quanta of these modes, i.e. phonons, as mediators 
of a coherent spin-spin interaction. However, in addition to 
the coherent dynamics, the phonons also provide an indirect 
coupling to the electromagnetic reservoir that leads to some 
collective dissipation on the spins. Our goal is to find suitable 
regimes where these collective effects allow for QC/QS. 

Tailoring the collective Liouvillian - To guide our search 
for the correct regime, note that the two-qubit gates for QC 
presently implemented in different laboratories [8] use nearly- 
resonant spin-dependent forces. These gates rely on the mo- 
tional coherence to suppress the residual spin-phonon entan- 
glement. Since the motional coherence is absent in our case, 
the only possibility would be to work in the far off-resonant 
regime for QC[16] or QS[24], where |^^| <C |4| < such 
that 5n = cOa — {cOq — COn) • We identify below the additional 
conditions to tailor the collective coherent/dissipative phonon- 
mediated processes in presence of sympathetic cooling . 

For the particular values considered below, the laser- 
cooling rate reaches values of Wn ~ 10~^a)„. In this case, 
the sympathetic cooling is very strong, and the vibrations 
reach the steady-state very fast. Hence, we can use a second 
Schrieffer- Wolff transformation to integrate out the phonons 
from Eq. (4), and obtain the following effective Liouvillian 

Act = ^eff(McT) = -i[i:^eff , Mcr] + %ff (l^a)^ (6) 
where {la = Trph{/i}. Here, the coherent Hamiltonian is 

= I (^'cT+ar + H.c.) + £ ifif af , 



which contains the phonon-mediated interactions of strength 
Jfj^, and an effective Zeeman spHtting BJ^. The coupHng to 
the phonons also introduces an indirect dissipation via 

hi 

+ £(n7 + r'^f)(2arAi^a; - a+ar^a - ^aofar), 



where Tf^^Tf^^ are the effective dissipation rates. To find the 
correct regime, we must compare the time-scales derived from 



where we have introduced 5„ = 5„ + ~ (^n )*}' ^« 

stands for the steady- state mean phonon number after the sym- 
pathetic cooling. From the above expression, it is clear that by 
tuning the ratio ^„ = W„ (/t„ + 1 ) / 1 5„ | , we can control whether 
the phonon-mediated interactions prevail over the dissipation 
^„ <C 1, or vice versa ^n^^- Below, we estimate the orders 
of magnitude for realistic values of an experimental setup, and 
describe how to use both regimes for QC/QS. 

Coherent and dissipative generation of entanglement - We 
consider the simplest scenario to test our scheme: a three-ion 
chain in a linear Paul trap (Fig. 1(b)). To use realistic param- 
eters, we take +Mg^^-+Mg^^-+Mg^^ as a prototype example, 
although we emphasize that the scheme also works for other 
species. We set the axial trap frequency to (0^/2% = 4.1MHz, 
which is possible by optimizing the trap voltages. The dipole- 
allowed transition \g) = \3Si/2) ^ k) = 13^/2) of ^Mg^"^, 
which is characterized by ^ 280. 3nm and r^/27r ^ 41.4 
MHz, shall be used for continuous sympathetic cooling. By 
applying an external magnetic field, we can encode the spins 
in a couple of Zeeman sublevels \F,M) of the ground- state 
manifold of +Mg2^ namely |t) = |2,2) and ||) = |3,3). This 
leads to a resonance frequency of (0^ /In ^ 1.79 GHz, and a 
negliglible decay rate of T(y/27t ^ 10~^^Hz. Finally, a pair 
of weaker laser beams drive the axial red- sideband transition 
through a far-detuned excited state in the 3/3/2 manifold of 
+Mg^^, such that the Lamb-Dicke parameter is T]^ ^0.16. 

The particular isoptopic mass ratio m^/nia ^ 0.96 implies 
that the axial vibrational modes are almost unchanged with 
respect to the homogeneous chain, a)„/27r^ {4.1,7.1, 10.1} 
MHz. In order to attain a wide range of possible values 
for the ratio we tune the Raman lasers closer to the so- 
called egyptian mode, such that the detunings are 5„/2;r G 
{6.2,3.2,0.3}MHz. Note that due to the large detuning from 
the remaining modes, the collective effects will be mediated 
by the egyptian mode. To sympathetically cool it, we place 
the cooling isotope at the middle of the chain, such that the 
standing-wave frequency is red-detuned (0^ = (oj + with 

= — r^/2. According to Eq. (5), this leads to a steady-state 
mean phonon number n^ =0.65 independent of the standing- 
wave Rabi frequency. Therefore, we can modify the Rabi fre- 
quency 0.2 < ^l^/rr < 2 in order to control the cooling rate 
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Figure 2. Coherent/dissipative generation of entanglement: (a) 

Coherent exchange dynamics (main panel) for the Liouvillian (4) 
specific to the egyptian mode (solid lines) and the effective descrip- 
tion (6) (circles and squares). Bell-state fidelity ^ = |(v^b|Mcj|V^b)| 
during a single SWAP (left lower panel), and optimization of the fi- 
delity for different gate times (right lower panel). Solid lines rep- 
resent the complete Liouvillian (4), whereas the stars follow from the 
effective description (6). (b) Dissipative dynamics (upper panel) un- 
der Eq. (4) (solid lines) and Eq. (6) (circles and squares). Bell- state 
fidelity ^ = | (0_ |/icj|0_) | at the steady state (lower panel). 



W3, and thus the ratio ^3, which shall allow us to explore 
regimes of either dominant collective dissipation or coherent 
interactions. We remark that the laser used for cooling ^Mg^^ 
will be highly detuned from the cooling transition of +Mg^^ 
(i.e. A/2;r ^ 2.7 THz), such that the induced decay rate for 
the considered regime lies in r^(^l^/A)^/ (2;r) < 10~^Hz, and 
will not affect the dynamics. This laser only contributes with 
off-resonant ac- Stark shifts that shall be considered later on. 
We now explore two possible applications for this setup. 

(a) Coherent generation of entanglement: Our goal is to use 
the coherent phonon-mediated interaction to generate entan- 
glement between the +Mg^^ ions. By setting = O.ISF^, 
we obtain a cooling rate of Ws/cos ^ 4.3 • 10~^, such that 
^3^1- 10~^, and the Hamiltonian part of the Liouvillian (6) 
thus dominates. Initializing the spin state in |i/Acy(0)) = Itiis). 
and setting Q^aVs ^ 10^3. such that the distance between 
the ions is an integer multiple of the effective Raman wave- 
length, we obtain the Bell state |i/Ab) = ^(Itils) - illits)) 
for ^ 4ms (Fig. 2(a)). We have considered a heating rate 
Fah ~ O.lphonon/ms for the rf trap. Even if the process is 
slower than the usual gates [8], it prevails over the phonon- 
mediated decoherence, and leads to fidelities above 98%. 

(b) Dissipative generation of entanglement: A different 
possibility is to exploit the collective dissipation to generate 
entanglement in the steady state. The idea is to set = 2F^, 
such that Ws/cos ^ 0.03 and ^3 ^ 2.3, and the dissipative part 
of the Liouvillian (6) becomes dominating. In this regime, we 
can exploit a superradiant/subradiant phenomenon [25]. Con- 
sidering the different decay channels = Vf^ ± F^f the 
symmetric and antisymmetric Bell states = ^^(Itils) ± 

lilts))' it is possible to set the ratio between the ions dis- 
tance and the Raman wavelength, such that either F^^ = 
or r^I^ = 0. We note that this phononic superradiance goes 
beyond the possibilities of the pioneer experiments of pho- 
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tonic superradiance with trapped ions [26]. In particular, start- 
ing from from 1 1/^(7(0)) = Itils) in the regime of T^I^ = 0, 
we can generate entangled states dissipatively. After set- 
ting Q^oTls ^ W3, we obtain a decoherence-free steady-state 
Bell state for r > ^ss ^ ^Ofis (Fig. 2(b)) with fidelities 
around 30%. We note that the ultimate limit of 50% cannot be 
achieved due to the thermal contribution to Eq. (6). However, 
schemes originally formulated for cavities [27], can also be 
applied in our trapped-ion setting to reach unit fidelities. 

Note that the typical time-scales of coherent and dissipative 
processes are 0.1- 10ms. In this regime, some other sources 
of noise may also become relevant. We show below that our 
protocol can be extended to make it robust against such noises. 

Other sources of noise - Fluctuating magnetic fields, laser 
intensities and thermal noise, lead to the dephasing term 

Hn=l^\{Bif + Fi{t))af. (7) 

i 

Here, BJ^ introduces noise via thermal fluctuations over the 
steady-state mean phonon number, and Fi{t) is a Markovian 
Gaussian random process that models the noise of external 
magnetic fields and uncompensated ac-Stark shifts, such as 
those provided by the cooling laser. This term leads to a de- 
phasing rate of Fd = ^ Jq dt' {Fj{t')Fj)st = 1/272, where we 
have assumed a local noise with a short correlation time Tc, 
and {Fj{t^)Fi)^t = ^e~^'/'^'^5;/. Let us note that that typical 
decoherence times are of the order T2 ~ 1-lOms (F^ ^ 0.05- 
0.5 kHz), and that for the regimes considered above, one ob- 
tains BJ^/Hn ~ 0.3-3 kHz. Comparing these values to the 
above time-scales (~0.1-10ms), the need for a scheme to 
actively decouple from these noise sources becomes rather 
clear. A partial solution are the so-called clock states, such as 
It) = |2, 1) and |;) = |3, 1) at a field of = 213G for ^Mg^\ 
which become insensitive to linear Zeeman shifts. However, 
we still need to mitigate the other sources of dephasing by a 
different mechanism. As a byproduct, this decoupling mech- 
anism yields also a tool to tailor the collective Liouvillian. 

The standard approach for prolonging the coherence of a 
system consists of a sequence of pulses that refocus the in- 
duced dephasing. This technique is known as pulsed dynami- 
cal decoupling [28], and has already been applied to ions [29]. 
Another approach known as continuous dynamical decou- 
pling, produces a similar effect by using continuous drivings, 
and has also been successfully applied to ions [30]. We apply 
a strong continuous driving resonant with the spins [31], such 
that the bare spin Hamiltonian (1) now reads 

Here, the driving frequency is set to C0(^ = cOa, and its Rabi 
frequency \Q.(\\^COa. Note that the considered resonance fre- 
quencies lie in the microwave regime (Oq/Iti^ 1 GHz, such 
that the driving can still be strong enough to fulfill 

{W„4}<|^^d|, Fd,T-i<|^ld|. (8) 

In this regime, we apply a modified Schrieffer- Wolff trans- 
formation to obtain an new master equation = ^eff(McT)' 
where the collective Liouvillian is modified as follows 

^eff(McT) = -i[^eff + ^n,McT]+^eff(McT) + ^n(McT). (9) 



The coherent Hamiltonian above describes a collection of 
spins interacting according to the so-called Ising model 

i>j in 

whereas the phonon-mediated dissipation corresponds to a 
collective dephasing in the eigenbasis of the strong driving 

ij 

Let us emphasize that the collective interaction strengths and 
dissipation rates, which are expressed as follows 

r 52 + ^2 2' ^ 52+^2 2^^^+2J^ 

lead to a very similar control parameter ^„ = (/t„ + ^ ) / 1 5„ | . 
Accordingly, under the same assumptions considered above, 
we can interpolate between regimes where the coherent Ising 
interactions dominate over the collective Ising dissipation, or 
viceversa. In addition to the new range of possibilities offered 
by this collective Liouvillian, note that the noise terms in (7) 
contribute to the new Liouvillian (9) as follows 

Hn=Y,^aaGf, #n(Ma)=I^fd(cjfAlaCjf-Ala), 

i i 

where Cl^ = {Bf^BnTc ^2r^)/2£l^T^, and fd = Fd/(Ild^c)^. 
According to the constraints (8), these noisy terms are clearly 
suppressed by a sufficiently strong driving. To be more spe- 
cific, for noise correlation times on the order of Tc = 10~^72 
and detunings 5„/2;r ~0. 1-1 MHz, it suffices to apply driving 
Rabi frequencies of Q^/27l ^ 10 MHz to reduce both sources 
of noise by more than two orders of magnitude. Therefore, 
we are decoupling from the noise very efficiently, while pre- 
serving the collective part of the Liouvillian for QC and QS. 
In contrast to Fig. 2(a), we remark that the new Liouvillian (9) 
allows for the coherent generation of all four Bell states with 
similar fidelities and gate times. Besides, as a consequence 
of our control over the collective dephasing, any of these Bell 
states is made decoherence-free by tuning the ratio of the ion 
distance and the laser wavelength appropriately. 

Many -body physics- The collective Liouvillians in Eqs. (6) 
and (9) define our toolbox for dissipation-assisted QS. So 
far, we have considered quantum information processing by 
means of an isolated vibrational mode that mediates the col- 
lective dynamics. In the many-ion scenario, this would lead to 
long-range interactions that define a fully-connected spin net- 
work. For surface traps (Fig. 1(a)), the presence of magnetic 
frustration already leads to rich phase diagrams [32]. How- 
ever, to achieve strongly-correlated models, it is better to work 
with a full vibrational branch of a small frequency width [24] . 
The small ratios ^ 10~^ obtained for trapping frequencies 
COt/271 ^ lOMHz, indicate that the coherent dynamics can be 
dominating also in this case, while minimizing the heating. 
This would allow for the QS of exotic models [33] with 3- 
body interactions or topological order in surface ion traps. 
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These conventional many-body QS focus on collective 
properties of the Hamiltonian while minimizing the influ- 
ence of the environment. However, the dissipative dynam- 
ics may also lead to interesting many-body phenomena [34] . 
The possibility to control the relative strength of the co- 
herent/dissipative parts in (6) and (9) is very appealing in 
this respect. We note that for intermediate driving strengths 
Q,(i, competing dissipative terms in (9) can lead to purely- 
dissipative quantum phase transitions. 

Conclusions and Outlook - We have proposed a scheme 
based on sympathetic cooling to overcome the anomalous 



heating in surface traps, while allowing for QC and QS. In 
addition, the sympathetic cooling becomes a tool to tailor 
the collective effects of the Liouvillian. By controlling a 
single parameter, namely the laser-cooling power, we have 
shown how the Liouvillian interpolates between regimes of 
dominating coherent interactions/collective dissipation, both 
of which allow for generation of entanglement. Moreover, we 
have pointed out how this control can be exploited for coher- 
ent/dissipative many -body QS. 
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